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Dark Matter from i?^-gravity 



Jose A. R. Cembranos 
William I. Fine Theoretical Physics Institute, University of Minnesota, Minneapolis, 55455, USA 

The modification of Einstein gravity at higli energies is mandatory from a quantum approach. In 
this work, we point out that this modification will necessarily introduce new degrees of freedom. 
We analyze the possibility that these new gravitational states can provide the main contribution to 
the non-baryonic dark matter of the Universe. Unfortunately, the right ultraviolet completion of 
gravity is still unresolved. For this reason, we will illustrate this idea with the simplest high energy 
modification of the Einstein-Hilbert action: _R^-gravity. 
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INTRODUCTION 

Different astrophysical observations agree that the 
main amount of the matter content of our Universe is 
in form of unknown particles that are not included in the 
Standard Model (SM). Typical candidates to account for 
the missing matter can be found in well motivated ex- 
tensions of the electro- weak sector. However, there is a 
fundamental sector in our model of particles and interac- 
tions, where the introduction of new degrees of freedom 
is not only well motivated, but absolutely necessary. 

The non-unitarity and non-renormalizability of the 
gravitational interaction described by the Einstein- 
Hilbert action (EHA) demands its modification at high 
energies. The main idea of this work is to realize that this 
correction cannot be accomplished without the introduc- 
tion of new states; these states will typically interact with 
SM fields through Planck scale suppressed couplings and 
potentially work as dark matter (DM). 



71 -GRAVITY 

In spite of many and continuous efforts, the ultraviolet 
(UV) completion of the gravitational interaction is still 
an open question. In these conditions, it is difficult to 
make general statements about its phenomenology. We 
will adopt a very conservative and minimal approach in 
order to capture the fundamental physics of the prob- 
lem. The simplest correction to the EHA at high energies 
is provided by the inclusion of four-derivative terms in 
the metric that preserve the general covariance principle. 
The most general four-derivative action supports, in ad- 
dition to the usual massless spin-two graviton, a massive 
spin-two and a massive scalar mode, with a total of eight 
degrees of freedom (in the physical or transverse gauge 
In fact, four-derivative gravity is renormalizable. 
However, the massive spin-two gravitons are ghost-like 
particles that generate new unitarity violations, break- 
ing of causality, and inadmissible instabilities [sj . 

In any case, in four dimensions, there is a non-trivial 
four-derivative extension of Einstein gravity that is free of 



ghosts and phenomenologically viable. It is the so called 
i?^-gravity since it is defined by the only addition of a 
term proportional to the square of the scalar curvature to 
the EHA. This term does not improve the UV problems 
of Einstein gravity but illustrates our idea in a minimal 
way. In fact, i?^-gravity only introduces one additional 
scalar degree of freedom, whose mass mo is given by the 
corresponding new constant in the action, as one can see 
in Eq. H]): 
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where Mpi = (SttG^v)"^/^ ^ 2.4 x 10^** GeV, A ~ 
2.3 X 10"'^ eV, and the dots refer to higher energy correc- 
tions that must be present in the model to complete the 
UV behaviour. In this work, we will show that just the 
Action ([T]) can explain the late time cosmology since the 
first term can account for the dark energy (DE) content, 
while the third term is able to explain the dark matter 
(DM) one. The first term is just the standard cosmolog- 
ical constant, that we will neglect along this work. We 
will focus on the new phenomenology that introduces the 
third term when it can be identified with the observed 
DM. i?^-gravity modifies Einstein's Equations (EEs) as 
d, 0| (following notation from @ ) : 
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where Iai3f_cu = (30/33^"^ - ffa^ff/Ji/)- The new terms do 
not modify the standard EEs at low energies except for 
the mentioned introduction of a new mode. In fact, if we 
impose to preserve standard gravity up to nuclear den- 
sities or Big Bang Nucleosynthesis (BBN) temperatures, 
the constraints on toq are just mo 10^^^ eV. In this 
case, the new terms are expected to be negligible at den- 
sities lower than p ~ Tbbn^ ~ (100 MeV)'' or curvatures 
lower than Hbbn^ Ibbn'^/Mpi^ - (lO^^^ eV)^. In 
this paper, we will discuss in detail the restrictions and 
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FIG. 1: Evolution of the scale factor of the Universe: a(t) 
as function of time t (time is normalized to the age of the 
Universe to ^ 4.3 x lO" s [24[, and a(to) = !)■ The standard 
evolution is modulated by a coherent oscillation. Although 
this oscillation has a very small amplitude, it has associated 
a high frequency (determined by the mass of the scalar mode, 
mo — 1 eV in the figure) and it can constitute the observed 
amount of DM (we do not take into account the inhomo- 
geneities coming from the clusterization process). 



possible signatures of the model. We will see that the 
ones already discussed are not dominant. 

It is straight forward to check that the metric g^i, = 
[1 + ci sin(r7ioi)]r7^i/ is solution of the hnearized Eq. ([2]), 
i.e. for ci ^ 1, without any kind of energy source. In 
this work we will argue that the energy stored in such 
oscillations behaves exactly as cold DM and can explain 
the missing matter problem of the Universe (see Fig. [T]). 
We want to emphasize that this new mode of the metric 
is an independent degree of freedom that eventually will 
cluster and generate a successful structure formation if it 
is produced in the proper amount. 



INTERACTION WITH STANDARD MODEL 
PARTICLES 



perturbatively as 
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where g^i, is its classical background solution, h^i, takes 
into account the standard two degrees of freedom associ- 
ated with the spin-two (traceless) graviton, and (j) corre- 
sponds to the new mode. This scalar field has associated 
a canonical kinetic term with the mass mo as we have 
already commented. 

We will deduce the couplings of this scalar graviton 
with the SM fields by supposing that gravity is minimally 
coupled to matter (in the Jordan frame). In such a case, 
(j) is linearly coupled to matter through the trace of the 
standard energy-momentum tensor: 



(4) 



It implies that the couplings with massive SM particles 
are given a tree level. In particular, the three body in- 
teractions are given by: 



trec-lcvcl 
(t>-SM 



1 



(5) 



with the Higgs boson ($), (Dirac) fcrmions (ip), and elec- 
troweak gauge bosons, respectively. In contrast with 
what has been claimed in previous studies, this field 
does couple to photons and gluons due to the confor- 
mal anomaly induced at one loop by charged fermions 
and gauge bosons. We find (following notation from [8]): 
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The particular value of the couplings {cem and cg) de- 
pends on the energy. We will be particularly interested 
on the coupling with photons, which leads to potential 
observational decays of </>. We will perform all the calcu- 
lations restricting ourselves to the content of the SM but 
the exact values of the couplings depend also on heavier 
particles, charged with respect to these gauge interac- 
tions, that may extend the SM at higher energies. 



In order to be quantitative, we need to write the action 
for the new scalar degree of freedom of the metric in a 
canonical way. This work can be done directly ,2] (what it 
is called inside the Jordan frame) or through a conformal 
transformation of the metric [3] (what it is known as the 
Einstein frame). As we will work in the limit in which 
R ^ TOq, in both cases, the metric can be expanded 



ABUNDANCE 

In principle, the above interactions of the scalar mode 
with the SM could produce a thermal abundance of (/> at a 
very early stage of the Universe. However, it is expected 
that higher order corrections to Action ([T]) will be im- 
portant at this point. In fact, it will typically take place 
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at temperatures T » Kq = ^/Mplmo, when we need to 
know the UV completion of the gravitational theory to 
study its dynamic. Nevertheless, there is at least another 
abundance source for this scalar mode that can be com- 
puted with Eq. (1). In the same case that other bosonic 
particles, such as axions this field may have associ- 
ated big abundances through the so called misalignment 
mechanism. There is no reason to expect that the ini- 
tial value of the scalar field ) should coincide with the 
minimum of its potential (0 = 0) if H(T) ^ toq- Below 
the temperature Ti for which 3H(Ti) ~ toq, ((> behaves 
as an standard scalar. It oscillates around the minimum. 
These oscillations correspond to a zero-momentum con- 
densate, whose initial number density: too'/'i/2 
(where (fii = ^ {(j){Ti)^) ), will evolve as the typical one 
associated to standard non-relativistic matter. 

Taking into account that the number density of scalar 
particles scales as the entropy density of radiation (s = 
27r^gsirf /45) in an adiabatic expansion, we can write: 
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where pcrit = 1.0540 x lO'^eVcm"'^ is the critical den- 
sity, So = 2970 cm^^ is the present entropy density of 
the radiation, and 7^1 is the factor that this entropy has 
increased in a comoving volume since the onset of scalar 
oscillations. 

If we supposed a radiation dominated universe at Ti 
{3Hi = 7r(gei/10)^/2TiVMpi), we can estimate Ti by 
solving mo = 3Hi{Ti): 



Ti ~ 15.5 TeV 
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where g^i (5si) are the effective energy (entropy) number 
of relativistic degrees of freedom at Ti . We see that initial 
values for the scalar field of order of 0i ~ lO^^ GeV 
can lead to the non-baryonic DM (NBDM) abundance 
depending on the rest of parameters and the early physics 
of the Universe (see Fig. [2]). We can check that this 
result is consistent with a perturbative treatment of the 
background metric | |A5p,y/5p,y| | < 10^^ for the entire 
computation. 



SIGNATURES AND CONSTRAINTS 



On the other hand, Eqs. (I4|5p imply that the new 
scalar graviton mediates a standard attractive Yukawa 
force between two non-relativistic particles of masses Ma 
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with a = 1/3 0. The non-observation of such a force by 
torsion-balance experiments requires 13, 11|: 



Too > 2.7 X 10"^eV at 95 % c.l. 



(11) 



This is the most constraining lower bound on the mass of 
the scalar mode and it is independent of its misalignment 
or any other supposition about its abundance. 

On the contrary, depending on its abundance, toq is 
constrained from above. It is particularly interesting the 
decay in e'^e^ since it is the most constraining if (j) con- 
stitutes the total NBDM. From (O, it is possible to cal- 
culate the (j) decay rate into a generic pair fermion anti- 
fermion. We find: 
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where is the number of colors, and = 1 or = 2 
depending if the fermion is Dirac or Majorana type. In 
particular, for the electron-positron decay {Nc = 1, /.t = 1 
and defining — TOo/(2TOe)): 
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Restrictions are set by the observations of the SPI spec- 
trometer on the INTEGRAL (International Gamma-ray 
Astrophysics Laboratory) satellite, which has measured 
a 511 keV line emission of 1.05 ± 0.06 x 10~^ photons 
cm~2 g-i from the Galactic center (GC) [l^, confirming 
previous measurements. This 511 keV line flux is fully 
consistent with an e+e^ annihilation spectrum although 
the source of the positrons is unknown. 

If Too > 1.2 MeV, the scalar field cannot constitute the 
total local DM since we should observe a bigger excess 
of the 511 line coming from the GC. On the other hand, 
decaying DM (DDM) has been already proposed in dif- 
ferent works as a possible source of the inferred positrons 
if its mass is lighter than Mnpivr IP.-^.^ Y Jl'^ 
decay rate in e+e" verifies [14l. Ill 
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Mddm 



[(0.2 - 4) X 10^7 s MeV] ^ . (14) 



The most important uncertainty for this interval comes 
from the dark halo profile, although a cuspy density is 
definitely needed (with a inner slope 7 > 1.5 [l^])- H 
Too is tunned to 2TOe with an accuracy of 5 — 10 %, 
the line could be explained by i?2-gravity. The same 
gravitational DM can explain the 511 line with a less 
tunned mass (up to Too ~ 10 MeV) if the misalignment 
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is ^ 10^ GeV, i.e. with a lower abundance (See Fig. 
[2]). If mo ^ 10 MeV, the gamma ray spectrum originated 
by inflight annihilation of the positrons with interstellar 
electrons is even more constraining than the 511 keV pho- 
tons [H. 

On the contrary, if mo < 2 mg , the only decay channel 
that may be observable is the decay in two photons. We 
find: 
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By taking into account all SM charged particles and as- 
suming (f) to be much lighter than all of them: cem = 
11/3, and 
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As it has been discussed in detail in [3, US]) if mo ^ 
1 MeV, it is difficult to detect these gravitational decays 
in the isotropic diffuse photon background (iDPB). The 
gamma-ray spectrum at high Galactic latitudes can have 
contributions from Galactic and extragalactic sources, 
but it seems well fitted at i?^ ^ 1 MeV by assuming Ac- 
tive Galactic Nuclei (AGN) as main sources. The spec- 
trum observed by COMPTEL (-the Compton Imaging 
Telescope- over the energy ranges 0.8 — 30 MeV [20||), 
SMM (-the Solar Maximum Mission- for 0.3-7 MeV [2l[) 
and INTEGRAL (5-100 keV faU hke a power law, 
with dN/dE ~ E''^ '^ and it will dominate any pos- 
sible signal from i?^-gravity if mo ^ 1 MeV. 

However, a most promising analysis is associated with 
the search of gamma-ray lines at E^ = 1^0/2 from local- 
ized sources, as the GC. The iDPB is continuum since 
it suffers the cosmological redshift. However, the mono- 
energetic photons originated by local sources may give 
a clear signal of i?^-gravity. INTEGRAL has performed 
a search for gamma-ray lines originated within 13° from 
the GC over the energy ranges 0.08 — 8 McV. It has not 
observed any line below 511 keV up to upper flux limits 
of 10~^-10~^ cm~^ s~^, depending on line width, energy. 



and exposure 23[. Unfortunately, these flux limits are, 
at least, one order of magnitude over the expected fluxes 
from (j) decays with mo ^ 1 MeV, even for cuspy ha- 
los. The photon flux originated by i?^-gravity depends 
on mo as ^Ey=mo/2 oc mQ. This strong dependence im- 
plies that only the heavier allowed region could be de- 
tected with reasonable improvements of present experi- 
ments fil. 



CONCLUSIONS 




mo(eV) 

FIG. 2: Parameter space of the model: mo is the mass of the 
new scalar mode and is its misalignment when 3H ~ mo 
(we assume gei ~ Qsi — 106.75, and 7si ~ 1). The left side 
is excluded by modifications of Newton's law. The right one 
is excluded by cosmic-ray observations. In the limit of this 
region, _R'^-gravity can account for the positron production 
in order to explain the 511 keV line coming from the GC 
confirmed by INTEGRAL [H (up to mo ~ 10 MeV). The 
upper area is ruled out by DM overproduction. The diagonal 
line corresponds to the NBDM abundance fitted with WMAP 
data fil. 



low energy phenomenology of the studied scalar mode is 
ubiquitous in high energy corrections of the EHA com- 
ing from string theory, supersymmetry or extra dimen- 
sions (reproduced in form of dilatons, radions, gravis- 
calars or other moduli fields) . The instability of the de- 
duced DM predicts a deviation from EEs at an energy 
scale: 1 TeV ^ Ac ^ 10^ TeV. Consequences for hier- 
archy interpretations, baryogenesis or inflation deserve 
further investigations. 
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In conclusion, we have studied the possibility that the 
DM origin resides in UV modifications of gravity. Al- 
though our results may seem particular of i?^-gravity, the 
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